function [adjMatrix,nodeNames] = getDAG()
nodeNames{1} = 'N0_7muVerMo';
nodeNames{2} = 'SubjVertMo';
nodeNames{3} = 'QGVertMotion';
nodeNames{4} = 'CombVerMo';
nodeNames{5} = 'AreaMeso_ALS';
nodeNames{6} = 'SatContMoist';
nodeNames{7} = 'RaoContMoist';
nodeNames{8} = 'CombMoisture';
nodeNames{9} = 'AreaMoDryAir';
nodeNames{10} = 'VISCloudCov';
nodeNames{11} = 'IRCloudCover';
nodeNames{12} = 'CombClouds';
nodeNames{13} = 'CldShadeOth';
nodeNames{14} = 'AMInstabMt';
nodeNames{15} = 'InsInMt';
nodeNames{16} = 'WndHodograph';
nodeNames{17} = 'OutflowFrMt';
nodeNames{18} = 'MorningBound';
nodeNames{19} = 'Boundaries';
nodeNames{20} = 'CldShadeConv';
nodeNames{21} = 'CompPlFcst';
nodeNames{22} = 'CapChange';
nodeNames{23} = 'LoLevMoistAd';
nodeNames{24} = 'InsChange';
nodeNames{25} = 'MountainFcst';
nodeNames{26} = 'Date';
nodeNames{27} = 'Scenario';
nodeNames{28} = 'ScenRelAMCIN';
nodeNames{29} = 'MorningCIN';
nodeNames{30} = 'AMCINInScen';
nodeNames{31} = 'CapInScen';
nodeNames{32} = 'ScenRelAMIns';
nodeNames{33} = 'LIfr12ZDENSd';
nodeNames{34} = 'AMDewptCalPl';
nodeNames{35} = 'AMInsWliScen';
nodeNames{36} = 'InsSclInScen';
nodeNames{37} = 'ScenRel3_4';
nodeNames{38} = 'LatestCIN';
nodeNames{39} = 'LLIW';
nodeNames{40} = 'CurPropConv';
nodeNames{41} = 'ScnRelPlFcst';
nodeNames{42} = 'PlainsFcst';
nodeNames{43} = 'N34StarFcst';
nodeNames{44} = 'R5Fcst';
nodeNames{45} = 'Dewpoints';
nodeNames{46} = 'LowLLapse';
nodeNames{47} = 'MeanRH';
nodeNames{48} = 'MidLLapse';
nodeNames{49} = 'MvmtFeatures';
nodeNames{50} = 'RHRatio';
nodeNames{51} = 'SfcWndShfDis';
nodeNames{52} = 'SynForcng';
nodeNames{53} = 'TempDis';
nodeNames{54} = 'WindAloft';
nodeNames{55} = 'WindFieldMt';
nodeNames{56} = 'WindFieldPln';
adjMatrix = zeros(length(nodeNames));
adjMatrix(27,49) = 1;
adjMatrix(27,48) = 1;
adjMatrix(24,36) = 1;
adjMatrix(35,36) = 1;
adjMatrix(22,31) = 1;
adjMatrix(30,31) = 1;
adjMatrix(21,22) = 1;
adjMatrix(5,13) = 1;
adjMatrix(9,13) = 1;
adjMatrix(12,13) = 1;
adjMatrix(28,30) = 1;
adjMatrix(29,30) = 1;
adjMatrix(15,20) = 1;
adjMatrix(16,20) = 1;
adjMatrix(5,9) = 1;
adjMatrix(8,9) = 1;
adjMatrix(27,45) = 1;
adjMatrix(31,42) = 1;
adjMatrix(36,42) = 1;
adjMatrix(40,42) = 1;
adjMatrix(41,42) = 1;
adjMatrix(21,24) = 1;
adjMatrix(23,24) = 1;
adjMatrix(27,51) = 1;
adjMatrix(5,21) = 1;
adjMatrix(13,21) = 1;
adjMatrix(19,21) = 1;
adjMatrix(20,21) = 1;
adjMatrix(27,50) = 1;
adjMatrix(27,32) = 1;
adjMatrix(27,28) = 1;
adjMatrix(16,19) = 1;
adjMatrix(17,19) = 1;
adjMatrix(18,19) = 1;
adjMatrix(13,15) = 1;
adjMatrix(14,15) = 1;
adjMatrix(4,5) = 1;
adjMatrix(27,56) = 1;
adjMatrix(38,40) = 1;
adjMatrix(39,40) = 1;
adjMatrix(27,53) = 1;
adjMatrix(27,52) = 1;
adjMatrix(27,47) = 1;
adjMatrix(27,46) = 1;
adjMatrix(6,8) = 1;
adjMatrix(7,8) = 1;
adjMatrix(26,27) = 1;
adjMatrix(37,43) = 1;
adjMatrix(42,43) = 1;
adjMatrix(27,55) = 1;
adjMatrix(27,41) = 1;
adjMatrix(32,35) = 1;
adjMatrix(33,35) = 1;
adjMatrix(34,35) = 1;
adjMatrix(15,25) = 1;
adjMatrix(15,17) = 1;
adjMatrix(16,17) = 1;
adjMatrix(27,54) = 1;
adjMatrix(25,44) = 1;
adjMatrix(43,44) = 1;
adjMatrix(27,37) = 1;
adjMatrix(10,12) = 1;
adjMatrix(11,12) = 1;
adjMatrix(1,4) = 1;
adjMatrix(2,4) = 1;
adjMatrix(3,4) = 1;
P = topologicalPermutation(adjMatrix);
adjMatrix = adjMatrix(P,P);
nodeNames = nodeNames(P);
